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ABSTRACT 


Many  monopulse  angle  tracking  radar  systems  are  made  up  using 
nonlinear  components,  such  as  limiters  and  logarithmic  amplifiers. 

This  report  describes  several  areas  of  work  that  all  relate  to  mono¬ 
pulse  angle  tracking  radars  having  at  least  one  limiter  in  the  signal 
processor.  In  one  area  the  detailed  angle  tracking  accuracy  was  de¬ 
termined  for  a  sum-and-difference  amplitude  monopulse  system  having 
a  sum  channel  limiter  prior  to  the  angle  error  detector.  The  analysis 
utilized  a  general  accuracy  expression  valid  for  an  arbitrary  angle 
tracking  system;  the  expression  was  developed  as  a  second  area  of 
effort.  A  third  area  of  work  centered  on  evaluating  coefficients  that 
enter  Into  the  accuracy  results  derived  in  the  first-described  area. 

A  fourth  work  area  Involved  computation  of  the  accuracy  of  a  Rhodes 
Class  III  monopulse  (with  two  limiters);  this  work  was  Incomplete  at 
the  project  end  and  remains  for  future  work. 
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1.0  INTRODUCTION 


Accuracy  of  angle  tracking  radar  systems  when  noise  is  present  has 
been  the  subject  of  many  analyses.  Until  the  earlier  phase  of  this 
study  [I]1"  most  work  did  not  employ  the  rigorous  analytical  methods 
currently  available.  In  [1]  rigorous  methods  were  used  to  develop  an 
accuracy  formula  for  the  mean-squared  angle  tracking  error  of  four  mon¬ 
opulse  systems.  These  four  systems  could  be  categorized  as  "linear"  be¬ 
cause  their  architectures  did  not  involve  nonlinear  signal  processing 
components,  such  as  limiters  or  logarithmic  amplifiers. 

It  is  the  purpose  of  this  report  to  describe  the  results  of  a  rig¬ 
orous  study  of  monopulse  angle  tracking  accuracy  applicable  to  systems 
with  nonlinear  elements.  The  reported  results  apply  only  to  architec¬ 
tures  having  a  limiter  or  limiters  present.  Architectures  having  log¬ 
arithmic  amplifies  (such  as  the  Rhodes  class  II  system  [2])  constitute 
an  area  of  possible  future  research. 

This  report  describes  work  that  falls  mainly  in  four  specific  areas. 
First,  the  tracking  accuracy  Is  found  for  a  sum-and-difference  channel 
form  of  amplitude  monopulse  system  when  the  sum  channel  contains  limit¬ 
ing;  this  is  probably  one  of  the  most  popular  architectures  in  use  to¬ 
day.  The  limiting  can  occur  because  the  sum  channel  input  to  the  angle 
error  (phase)  detector  is  saturated,  or  because  a  separate  limiter  is 
used  prior  to  the  angle  detector.  In  either  case  the  effect  can  be 
modeled  as  a  limiter  prececding  the  angle  error  detector.  This  type  of 
system,  which  will  be  called  system  5,  is  analyzed  in  Section  3.0. 

The  second  area  of  analysis  concerns  a  general  expression  for 

X 

'  References  are  cited  by  bracketed  number  and  are  listed  at  the  end 
of  the  report. 


the  angle  accuracy  of  an  arbitrary  system,  which  Includes  those  hav¬ 
ing  nonlinear  elements.  The  resulting  expression  formed  the  basis  for 
the  accuracy  results  derived  for  system  5. 

The  third  area  of  work  involved  the  numerical  evaluation  of  some 
coefficients  that  appear  in  the  accuracy  expressions  of  system  5.  This 
work  was  only  partially  completed  at  the  end  of  the  study  because  of 
some  computer  programming  difficulties. 

The  fourth  major  area  of  effort  centered  around  the  accuracy  analy¬ 
sis  of  a  Rhodes  Class  III  monopulse  angle  tracking  system  (one  with  two 
limiters).  This  work  was  partially  completed  and  is  expected  to  be 
completed  through  future  effort. 

Finally,  a  number  of  other  minor  tasks  were  undertaken,  mostly  to 
support  the  above  mentioned  four  efforts.  The  most  significant  other 
task  involved  extending  the  general  analysis  of  Kelly  and  Hariharan  [3] 
for  the  general  response  of  an  ideal  limiter  when  its  input  is  the  sum 
of  a  bandpass  periodic  signal  and  bandpass  Gaussian  noise. 
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2.0  GENERAL  THEORY  FOR  ACCURACY 
fiF  ANGLE  TRACKING  RADAR 

2.1  Introduction 

In  all  angle  tracking  radars  there  Is  some  point  where  an  "error" 
voltage  is  de, eloped  that  is  used  to  force  the  antenna  patterns  to 
follow,  or  track,  the  target.  There  is  one  such  error  voltage  de¬ 
veloped  for  each  angle  coordinate.  The  error  voltages  may  be  the  re¬ 
sult  of  a  "linear"  system,  such  as  a  sum-and-difference  monopulse  pro¬ 
cessor  having  a  "linear"  error  detector  or  they  can  result  from  a  non¬ 
linear  system.  A  commonly  used  nonlinear  system  is  the  ordinary  sum- 
and-difference  processor  where  the  error  detector  (a  phase  detector)  is 
operated  with  one  input  (the  sum  channel  input)  saturated.  Another  ex¬ 
ample  would  be  a  Rhodes'  class  III  [2]  monopulse  receiver  having  a  limiter 
in  each  channel. 

Regardless  of  how  the  error  voltage  is  developed  it  contains  all 
the  information  available  to  the  tracking  system  about  the  target's 
location.  Since  it  also  contains  noise,  its  use  in  estimating  (tracking) 
target  angle  will  have  errors.  In  the  following  section  we  develop  a 
general  equation  for  the  mean-squared  angle  tracking  error  (variance 
for  a  zero-bias  system)  cf  a  monopulse  radar  that  is  valid  even  if  the 
system  contains  nonlinear  elements  in  its  signal  processor.^ 

2.2  Accuracy  Expression  j 

In  Figure  2-1  let  v  represent  the  output  of  some  monopulse  angle 
tracking  system  that  is  to  be  used  by  an  angle  servo  to  track  a  target 
in  an  angle  coordinate  denoted  by  x.  Because  angle  dynamics  pre  usually 
small  in  most  systems,  there  is  usually  some  form  of  lowpass  filtering 
(LPF)  of  v  prior  to  its  use.  Let  u  represent  the  filter's  output.  We 
shall  a-.sume  the  filter's  bandwidth  to  be  small  relative  to  the  system's 
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Fig.  2-1.  Arbitrary  system  generating  an  angle 
error  voltage. 


pulse  rate  (PRF)  but  otherwise  can  be  as  large  as  necessary. 


The  angle  tracker  can  operate  only  on  the  "signal"  part  of  u. 
Because  noise  is  present  we  could  take  its  effect  into  account  by  de¬ 
fining  the  signal  as  the  average  value  of  u  E(U).  [Here  E( * )  denotes 
the  statistical  average  taken  over  the  noise  modeled  as  a  sample  func¬ 
tion  of  a  random  process  U(t)].  This  definition  is  only  partially  ac¬ 
ceptable  because  E(U)  is  still  a  function  of  time.  It  is  therefore  the 
time  average  of  E(U),  denoted  by  t^{E(ll)},  that  defines  the  "signal" 
component  of  u.  Thus, 

exM{E(U)}  (2.2-1) 

nX  “  u  -  ex  -  u  -  \fi  (E(U)} .  (2.2-2) 


If  the  noise  nx ( t )  is  modeled  as  a  sample  function  of  a  random 
process  N(t),  the  time  averaged  autocorrelation  function  of  the  noise 
process  will  reduce  to 

Rn(t)  =  cfl(RN(t,  t  +  t)}  =  i/UE  [N(t)  N(t  +  x)]> 

*  (/HR^t,  t  +  t) }  -«fl{E[U(t)]}  JUE[U(t  +  t]}.  (2.2-3) 


Next,  let  x  represent  the  coordinate  of  possible  target  angles 

(relative  to  boresight)  and  let  xQ  be  a  specific  value.  For  values  of 

x  near  x„ 
o 


e„(x) «  ev(x  )  + 


dev 


x'  o' 


dx 


x=x. 


(x  -  x0). 


(2.2-4) 


An  error-free  measurement  of  u  would  equal  ex(x0).  With  noise  we  have 
an  arbitrary  u  and  e  (x  )  is  replaced  by  e  (x)  of  (2.2-4)  where  (x-x  ) 

a  O  X  0 

is  interpreted  as  the  measurement  error  in  xQ  that  is  associated  with 
the  noisy  value  of  u.  Hence, 
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(2.2-5) 


(2.2-6) 


x=x_ 


The  mean-squared  error  in  angle  is  defined  as 
o2[x)  A  t/? (E  [(x-xQ)2]}. 

From  (2.2-6)  using  (2.2-7)  and  (2.2-3)  with  t30: 
2,X  ufl  {E[x2(t)]}  <A{Rx(t,t)} 

(j  (x)a  — - - 


(2.2-7) 


r  de 

X 

1 

2 

rdex 

1 

dx 

O 

X 

II 

X 

, 

dx 

X 

II 

X 

o 

\ _ 

(2.2-8) 


or 


c/Uyt.t)}  -  (^(E[U(t)]}) 
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All  . 

kr‘.' 

’dex 

1  —  . 

dx 
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x»x„ 

0  «■* 

(2.2-9) 


Usually  we  are  most  interested  in  accuracy  near  boresight  where 
x  =  xQ  3  0. 

Other  forms  can  be  obtained  for  (2.2-9)  which  is  a  basic  result 
applicable  to  nearly  any  system.  If  we  let  *^(ui)  and  *^v(w)  denote  the 
Fourier  transforms  of  Ry(t)  and  Rv(t),  the  time  averaged  autocorrela¬ 
tion  of  v(t),  then  it  can  be  shown  [4],  p.  450,  that 

*?„(«)  *  I  HM  |2  (2.2-10) 

where  H(w)  is  the  transfer  function  of  the  LPF.  Thus, 
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3,0  ANGLE  ACCURACY  OF  A  MONOPULSE 
TRACKING  RADAR  HAVING  A 
SUM-CHANNEL  LIMITER 

3.1  INTRODUCTION 

In  this  section  we  shall  find  the  angle  tracking  accuracy  of  a  sum- 
and-difference  form  of  amplitude  comparison  monopulse  radar  when  the 
sum  channel  contains  a  limiter.  Accuracy  will  be  determined  through 
the  mean-squared  angle  error  (same  as  error  variance  when  the  system 
has  no  bias,  as  is  the  case  here).  Thus,  we  shall  basically  seek  the 
solution  of  (2.2-14)  for  the  system  of  interest. 

Figure  3-1  illustrates  the  block  diagram  of  the  system,  which  we 
call  System  5  (four  systems  were  previously  studied  in  [I]).  Some 
functions  have  been  omitted  from  Fig.  3-1  because  they  do  not  affect 
the  accuracy  results.  For  example,  the  down-converting  mixers  usually 
present  in  both  channels  only  affect  the  frequency  at  which  the  compon¬ 
ents  operate  but  do  not  affect  signal  forms  (shape,  bandwidth,  etc.)  so 
do  not  affect  accuracy.  In  a  similar  manner  the  automatic  gain  control 
(AGC)  loop  usually  has  a  range-gated  envelope  detector;  we  shall  assume 
the  AGC  loop  works  well  enough  that  it  does  not  affect  accuracy  and 
these  gate/ detector  functions  can  be  omitted. 

In  Fig.  3-1  we  make  several  assumptions.  The  noises  that  are 
present  at  the  inputs  to  the  matched  filters  (MF’s)  are  assumed  to  be 
white,  stationary,  zero-mean  and  Gaussian  with  equal  power  density 
spectra  c^/2.  The  amplifiers  (AMP's)  are  assumed  to  be  identical  in 
gain,  are  affected  by  the  AGC  loop  in  the  same  manner  and  any  band¬ 
shaping  they  might  contribute  to  the  channels  in  practice  is  accounted 
for  in  the  transfer  function  of  the  matched  filters.  Thus,  the  noises 
and  n^  ate  Gaussian,  bandpass,  and  zero-mean  with  identical  power 
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spectra.  Their  powers  are  also  equal  and  denoted  by  a  . 

The  above  assumptions  are  identical  to  those  made  in  an  earlier 
study  [1]  for  System  1.  In  fact,  both  siqnals  and  noises  are  identi¬ 
cal  up  to  the  outputs  of  the  MF’s.  From  [1]  we  have 


e2^  "  c (0 )  n=-« 


c(t-xT-n  Tr)  cos  [wQ(t  -  xT)  +  Qq]  (3.1-1) 


V  k  x  °° 

e4(t)  -  -  "  5,  c(t-tT-nTR)  cos  [»0(t  -  tt)  +  0Q]  (3J‘2) 

4  c(0)  o3  n 


C ( t )  =  P(S)  p(C  +  t)  dF. 


(3.1-3) 


Here  is  the  peak  amplitude  of  the  sum  channel  pulses  at  the  input 
to  the  limiter  (LIM);  the  shape  of  these  pulses  is  determined  by  c(t) 
which  is  related  to  the  shape  of  transmitted  pulses  determined  by  p(t). 
The  parameter  is  the  two-way  target  range  delay,  TR  is  the  system 
repetition  period  (pulse  repetition  rate  is  fD  =  1/TD),  w  is  the  radar's 
carrier  angular  frequency,  0Q  is  an  arbitrary  phase  angle,  x  is  the 
target's  error  angle  from  the  boresight  axis  in  the  coordinate  of  inter- 


,  A  °3 

k  =  T~ 
m  Gr 


DG/.v’x) 


(3.1-4) 


©3  is  the  beamwidth  of  the  one-way  antenna  sum  pattern  between  -3dB 
points  in  the  x  coordinate  direction,  GA(x)  is  the  one-way  x-coordinate 
antenna  voltage  difference  pattern,  and  G„  is  the  one-way  voltage  gain 
of  t.  sum  pattern  on  the  boresight  axis. 

*  Accuracy  is  developed  only  for  the  x  coordinate.  An  analysis  for  the 
space-orthogonal  coordinate  would  be  identical  in  form.  Numerical  re¬ 
sults  would  be  different  in  the  two  coordinates  only  if  parameters  (0., 
km,  etc.)  are  different  in  the  two  coordinate  directions.  J 


Noises  n2  and  n4  have  equal  power  spectra,  as  given  In  [1],  that 
can  be  put  In  the  form 

2  2  2 

(“)  (“)  ’  2c[o7  K-o>l  +  lp(u  +  “0)l  ]’  (3.1-5) 


2 

where 


c(t)  — ►|P(cj)j  ; 


(3.1-6) 


P(“)  denotes  the  Fourier  transform  of  p(t),  and  the  double-ended  arrow 

2 

represents  a  Fourier  transform  pair  [transform  of  c(t)  is  |P(w)|  ]. 

It  Is  clear  from  (3.1-1)  that  the  single-pulse  peak  signal  power 
to  average  noise  power  ratio  of  the  waveform  at  the  input  to  the  limit¬ 
er,  denoted  by  (S^/N),  is 

( S] /N )  -  V*  /2a2.  (3.1-7) 

The  relationships  derived  or  stated  in  this  section  form  the  link 
with  the  earlier  study  and  define  most  of  the  quantities  common  to 
both  studies.  To  arrive  at  the  desired  system  accuracy  we  must  now  com 
pute  the  mean  value  and  autocorrelation  function  of  the  system's  re¬ 
sponse  [v(t)  in  Figure  3-1]. 

3.2  Mean  Value  of  System  Response 
Frcr  :1gure  3-1  we  have 

v(t)  *  g(t)  z(t)  =  g(t)  w(t)  y(t).  (3.2-1) 

Since  w(t)  *  e^(t)  +  n^(t)  and  n^( t)  is  statistically  independent  of 
n2(t),  and  therefore  independent  of  y(t),  the  mean  value  of  the  process 
V(t),  having  v(t)  as  a  sample  function,  is 

E[V(t)]  *  g(t)  E[W(t)]E[Y(t)]  =  g(t)  e.(t)  E  [Y(t)].  (3.2-2) 


To  evaluate  E[Y(t)]  we  note  that  the  limiter  Is  described  by  Its 
transfer  characteristic  g(X)  where 


Y  *  g(X) 


(3.2-3) 


Thus,  If  f^(x)  ^  the  probability  density  function  of  the  limiter's 
input  process  X (t )  and  f^  (n2)  is  the  density  of  the  process  N2 ( t ) 


representing  the  noise  n2 ( t )  we  have 

WLm  fx<x) dx'[  n  f| 


N2(x  -  e2>  dx 


fN2(x  ■  e2>  dx  -  J_  n2 


f  (x  -  «2)  dx. 


(3.2-4) 


since 


fxu)  -  (X  -  e2). 


(3.2-5) 


By  substituting  for  fN  (•)  which  Is  a  zero-mean  Gaussian  density  with 
2  ^ 

variance  a  ,  and  using  simple  variable  transformations  (3.2-5)  will 


reduce  to 


E[Y]  ■  erf 


(— ) 

\/T  al 


where  erf  (•)  is  the  error  function  defined  by 

2  a  2 

erf(a)  «  e  ^  <*£. 


(3.2-6) 


(3.2-7) 


It  is  also  known  that 


erf(x) 


2  l 


k  2k  +  1 
(-1)  x 


3  A-  k=0  k!  (2k  +  1) 


(3.2-8) 


1 


By  using  (3.2-8)  and  (3.2-6)  in  (3.2-2)  we  have 


E[V(t)]  -  g(t)  e.(t) 


k*0  k! (2k+l ) 


2 

/T  a 


.  2k+l 

7  1 


(3.2-9) 


3.3  Time-Averaged  Mean  Value 
By  defining 


Per(t)“  s  c(t  -  tt  -  nTD) 

n«-« 


9  4  «0(t  -  tT)  -  n  Tp 


(3.3-1) 


(3.3-2) 


and  substituting  these  as  well  as  (3.1-7),  (3.1-1)  and  (3.1-2)  into 
(3.2-9),  the  time  average  of  (3.2-9)  can  be  written 

2/r  k,  *  °  a,  \  *  r  (-nk  1 


,  2k+2  v 

•<fl{ 9(t)  [-IfSPj  cos-f,)}. 


(3.3-3) 


Because 


cos-2(e) 


{  fe2)  *  2  ^(2km+2)cosKm2-2m)6]} 


(3.3-4) 


the  typical  term  In  (3.3-2)  is  the  time  average  of  the  product  of  a 
lowpass  function  and  a  bandpass  function.  Such  terms  will  be  nearly 
zero  because  cos  [2(k+l-m)9j  cycles  rapidly  during  the  times  the  low- 
pass  terns  exist.  Hence,  only  the  constant  term  in  (3.3-4)  gives  rise 
to  a  nonzero  result  in  (3.3-3).  The  remaining  time  averages  are  of  a 
periodic  function  (period  =  TR).  With  the  definitions 
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,v2 

A 

dk+l  3 

-V2 

2k+2 

B9-]  rect(x9) 

dt 

-(yM 

'i.  |t|  <  V2 

0,  IH  >  Tg/2 

k 

A  (-1) 

0 

(2k>!  /3k+l  ' 

r . 

‘ck  (k!)2  (k+1)!  22K  \  ai  / 

(3.3-5) 


(3.3-6) 


(3.3-7) 


we  can  finally  write  (3.3-3)  as 

«» ■  u  ■«  (+-)' 


(3.3-8) 


Next,  if  we  recognize 

•  -r1 - 


W  T 
»  L  R 


(3.3-9) 


2km  (S,/N) 


as  the  smallest  angle  error  variance  possible  for  any  monopulse  system 
[1],  the  derivative  required  in  (2.2-14)  is 

k 


dv^(Ervl> 

dx 


VE  dt 


x* 


0  [2"S,!„(x)  T«  ] 


TTT 


k=0 


fsi  r 

*ck  \~Wl  * 


(3.3-10) 


As  an  intermediate  point  in  computing  overall  system  accuracy  we 
use  (3.3-10)  in  (2.2-14)  to  write 

a  ?  (x)  - y  R  ?  [  iJHRv(t,  t+r)}  d  t 

o  (x)=  - 


ii'-fwr 


(3.3-11) 


3.4  Autocorrelation  of  System  Response 

To  complete  the  solution  of  (3.3-11)  we  must  find  Ry(t,  t+x). 

From  (3.2-1) 

Ry(t,t+x)  =  E  [V(t)  V(t+x)]  =  g(t)g(t+x)  E[W(t)W(t+x)  Y(t)Y(t+x)] 

a  g(t)g(t+T)  E[W(t)W(t+x)]  Ry(t,  t+x).  (3.4-1) 

Since  W(t)  =  e^(t)  +  N^(t)  and  the  noise  is  zero  mean,  (3.4-1)  reduces 
to 

Ry(t,t+x)  =  g(t)g(t+x)  R,^(x)  Ry(t,t+x)  (3.4-2) 

where  RN  (x)  is  the  autocori elation  function  of  the  stationary  noise 
n4 

process  N4(t).  In  writing  (3.4-2)  we  have  used  the  approximation 

e4(t)«0  for  nea --bores ight  tracking. 

By  inverse  transformation  of  (3.1-5),  RM  (x)  will  be 

4 

RN4(t)  *  °2  "flf}  COs(wOT)  “  °2  ~f{f}  cos  (V  (3.4-3) 

where  we  define 

0  £  t.  (3.4-4) 

Thus, 

Ry(t,t+x)  *  g(t)g(t+x)  a2  fjfj-  cos  (e, )  Ry(t,t+x).  (3.4-5) 

Our  main  remaining  task  is  to  solve  for  the  autocorrelation  function 

Ry ( * ,  * )• 

If  we  define  the  joint  probability  density  function  of  the  pro¬ 
cess  X(t)  at  times  t  and  t  +  x  as  f^Xj.Xg),  where  X1  -  X(t)  and 

*  X(f+x),  then  the  responses  are  denoted  Y1  =  Y(t)  *nd  Y 2  =  Y(t+x)  » 


V*-**)  *  Etv{t)yct+x):  .  E[yty2]  .  E[j  x-t  -|xUtT)|j 


-x,  X. 


A1 

TOT 


rx7  T7T  fx  (xl,x2)  dxl  dx2‘ 

|1|  I x2 1 


(3.4-6 


Next,  we  note  that  in  the  x^  x2  plane 


X1  x2 


1 ,  quadrants  1  and  3 


TX]T"  lx2i  }  “1  *  Quadrants  2  and  4. 


Furthermore, 


fx(xrx2)  =  fN2(xi  '  e]»  x2  ’  e2^ 


(3.4- 


(3.4-8 


A  A 

where  e^  *  e2(t)  and  e2  =  e2(t+r)  and  (*,*)  is  the  joint  density 
of  noise  process  N2(t)  at  times  t  and  t  +  t.  Because  these  noises 
are  Gaussian 


(a  )  *  2  /  ■  ■  , 
N2  2tto 


(  o'  -  2pa8  +  3^  ) 

exp  i '  17TT7)  i 


(3.4-9 


where  p  is  the  correlation  coefficient  given  by 


o  Rn  (t). 

& 


(3.4- 


By  using  (3.4-7)  through  (3.4-10)  in  (3.4-6)  it  can  be  shown 


RyU.t+r)  =  -1  +  2  [P(^)  -  Q(-gr)]  +  41  [^-,  5-,p],  (3.4- 

where  [5,  p.  936] 

k  -£^/2 

P(k)  =  — - —  f  e  d£  *  t  +  i  erf  (k//T  )  (3.4- 


We  separately  expand  the  terms  in  (3.4-16).  From  use  of  (3.2-8) 


we  expand  the  first  term  as 
e 


erf(!l_)erfp— U  L.  z 

l  /r  a)  \f2o  )  *  1=1 

i+k 

(-1)  e2 


CO 

z 

0  k=0 


i+k  2i  +1  %  2k+l,..  . 

(t)  e,  lt+T> 


(3.4-19] 


(2i+l)  i!  (2k+l)  k!  (2a2)  1+k+1 


If  (3.4-17)  and  (3.4-18)  are  used  in  the  second  term  of 
(3.4-16)  we  have 


>  pnt' 
4  E  p 


— ■U^MTTrW-r-M  4) 


n=0  (n+D 

-  Cn]/2  [nJ/2 


co  m  g  jjn+l(r)[e2(t)] 

Z  Z  nm&  r _ 


n-2m+2r 


n=0  m=0  4*0 

n-2£+2P 


r=0  p=0  n-m-£  +  r  +  P 


•  [e2(t+r)j 


(3.4-20) 


where 


m+£+r+p  - 

(-1/2)  (n! ) 


nmJ’  2tt  (n+1 )  I  m!  (n-2m)!  £!  (n-2£)!  r!  p! 


(3.4-21) 


Next,  we  use  (3.3-1)  and  (3.3-2)  to  write  (3.1-1)  as 
V„ 


e2(t)  =  —  Per(t)  cos  (9). 


(3.4-22) 


Therefore, 

V 

e2(t+r)  ■  Per (t+r)  cos  (9  +  6j)  (3.4-23) 

where  (3.4-4)  has  been  substituted.  After  using  these  two  expressions 
with  (3.4-20),  (3.4-19)  and  (3.1-7),  we  can  write  (3.4-16)  as 


uu  w 


E  . L-.VJL 


o  k«o  (21+D  i!  (2k+l)kl 
21+1  2k+l  21+1  2k+l 

[hpFJ 


»  [n]/2  [n]/2  *  »  /.  S,  \ 

*  4  I.  *  E  *  1  Enml  (2  ”5/ 


n-m-£+r+p 


n=0  m=0  £=0  r=0  p=0 

^+1  n‘2m  +  2 r 


TcCcif  f  Per(t)  *|  fPer(t+T)  1 

Lc(0)J  LTIorJ  L  c (0)  J 


n-2£  +  2p 


n+1  n-2m+2r  r  ^n-2£  +  2p 

[cos  (6^]  [cos  (e)]  [coste  +  e^J  *(3.4-24) 


3.5  Time-Averaged  Autocorrelation 

Our  principal  goal  is  to  solve  (3.3-11)  which  requires  that  we 
find  (/) (Ry(t,  t+r)}  by  tine-averaging  (3.4-5)  using  (3.4-24).  We  pro¬ 
ceed  first  with  the  double-sum  term  in  (3.4-24).  First  we  substitute 
the  term  into  (3.4-5)  ar.d  form  the  time  average.  By  grouping  terms 

involving  time  t  we  observe  that  cos^+^(0)  cos^  ^(0  +  0^)  is  periodic 

with  period  T^,  and  cycles  rapidly  with  t  [recall  that  0  is  a  function 

of  t  through  (3.3-2)]  while  other  factors  vary  slowly  (are  baseband). 

21+1 -  2k+l 

By  Fourier  series  expansion  of  cos  (0)  •  cos  (0+9^ )  only  the 

constant  term  (in  t)  leads  to  a  baseband  result  which  is  unfiltered  by 

the  output  lowpass  filter.  The  time  average  then  equals  this  constant 

term  times  the  time  average  of  other  baseband  time  terms.  The  overall 

result  is  a  double  sum  of  terms  that  are  periodic  in  the  variable  0j. 

Again,  a  Fourier  series  expansion  is  used  and  only  the  constant 


(baseband)  terms  are  preserved.  The  final  result  becomes 


where 


=  I 
1=0 


(3.5-1) 


and  Tg  is  the  duration  of  a  typical  gate  pulse. 

The  procedures  described  above  are  followed  to  develop  the  second, 
or  remaining,  part  of  (3.4-5)  needed  in  (3.3-11).  We  shall  omit  the 
details;  the  final  result  is 


(3.5-4) 
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00  [n]/2  [n]/2  00  =°  5  /S1 

Z  Z  Z  Z  Z  lWrp  Bnmjirp\-j^ 

A  n=0  m=0  1=0  r=0  p=o 

n-m-2+rt-p 

B  CO  CO  - 

I  Z  Hik  RAik 

i=0  k=Q 

'S1 
v  N  } 

i+k 

# 

(3.6-3) 

The  main  result  is  (3.6-1).  This  result  has  exactly  the  same 
form  as  the  tracking  error  variance  for  other  monopulse  [1]  and 
conopulse  [6]  radars.  Therefore,  the  presence  of  the  sum-channel 
limiter  has  not  altered  the  form  of  the  accuracy  formula  but  the  co¬ 
efficients  R^  and  Rg  have  become  much  more  complicated.  In  fact,  they 
are  now  functions  of  signal-to-noise  ratio  in  complicated  series  form. 


3.7  Coefficient  Evaluation 

In  many  cases,  as  in  this  study,  it  is  convenient  to  determine 
the  integrals  involved  in  the  expressions  (3.5-5)  and  (3.5-2)  for  co¬ 
efficients  Knm£v.p  and  Hik,  respectively,  by  numerical  integration. 
However,  it  is  possible  to  develop  series  solutions  for  these  co¬ 
efficients,  as  is  done  in  this  section. 

We  begin  by  noting  that  both  coefficients  involve  an  integral  of 
the  form 


n  I  K 

cos  (0)  cos  (0  +  0^ )d0 


(3.7-1) 


-TT 


where  I+K  =  even  integer.  We  first  use  trigonometric  identities  and 
the  binomial  series  as  follows 


r  t  k 

Q  ~  j  cos  (0)  [cos (0^ )  cos  (0)  -  sin^)  sin(©)]  d  e 

-7 T 


I+K-r  r 

cos  (0)  sin  (0)  d0 


(3.7-2) 


where 


a  «  cos (6j )  (3.7-3) 

b  *  s1n(0^ ) .  (3.7-4) 

T  x  |C  m  m 

Since  cos  (0)  Is  even  for  all  I+K-r  while  sin  (0)  is  odd  for  r  odd 


the  Integral  Is  zero  for  r  odd.  The  integral  is  therefore  nonzero 
for  only  even  values  of  r.  Let  [K]  *  largest  even  integer  not  exceed¬ 


ing  K  and  define 


r  -  2k  ,  0  *  k  -  [K]/2. 


(3.7-5) 


Then, 


[K]/2/ 


K'2k,  u\2k  f  » 


I+K_2k 

cos  (0)  sin'<K(0)d0.  (3.7-6) 


Because 


sin2k  (0)  -  [l-cos2(0)]k  -  E  (Jj)  [-cos2(0)]k  " 

n*0 


(3.7-6)  becomes 

ClO/2,  . 


(3.7-7) 


2k  *  /l,  \  k-n 

4  <-b)  I  ft)  H)  (6)1 


I+K-2n 


*»  (3.7-8) 

The  power  I+K  -  2n  is  even  for  all  I,  K,  and  n  so  a  series  expansion 
for  [cos (9)3  k  2,1  can  be  substituted  [7,  p.  25]  to  evaluate  the  integ¬ 
ral  easily.  We  h=>ve 


'  2T  {(%")*  2  (?)  <=«  [2(m-q>0]J 


(3.7-9) 
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so  (3.7-8)  becomes 


Q  -  f'2 «)  »K-2k  >2k  j0  (nk)  (-Dk-n_^ _ (\:i :  y 

k=0  'tk/  n~°  .I+K  -  2n  \ - 2 - 7 


(3.7-10) 


Q  -  2n 


[K]/2  k 
Z  E 
k=0  n= 


£![)(£)  (-i)k'nc«K'Zk(0,)  sin2k(e1)^l*|Ja) eI+k‘‘ 


(3.7-11) 


We  apply  (3.7-11)  to  the  solution  cf  (3.5-6).  After  noting  that 


in  this  case 


I  »  n  -  2m  +  2r 
K  *  n  -  2i  +  2p, 


(3.7-12) 

(3.7-13) 


we  substitute  (3.7-11)  into  (3.5-6)  which  will  then  contain  series  terms 
involving  Integrals.  The  integrand  of  a  typical  term  is 

[cos  (6^ )]2^n -£+P+l-k)  .[sfn(0j )]Zk.  The  cosine  factor  is  even  and 
its  power  Is  2  or  more.  The  sine  factor  is  expanded  according  to 

[sfn(0,)]2k  -  [1  -  cos2(e1)]k  ■  i  ^  (-1)1  (J)  cosZ1(e,)  (3'7''4) 

and  the  resulting  Integrand  is  expanded  using  known  series 

[7,  p.25,  eq.  1 .320(5. )]  and  integrated  easily.  The  final  result  is 


[n-2M2p)/2  k  k.s 

B_,„,  *  £  1  H>  \ 

nmj,rp  (<=  q  s=  0 


*  .k-s  fn-2£+2p)/ k y2(n-m-£+r+p-s)) 
£  M)  \  2k  AsA  n-m-Jl+r+p-s  / 


1  E  ,  ,  J  A)  J_(2J\ 

^n-m-£+r+p-s  ^_0  "  1  4J  ^ '  » 


(3.7-15) 


I 


/ 


J  =  n-  £  +  p  +  l  -  k  +  1 


(3.7-1 


u 


Coefficients  Knmirp  result  readily  by  using  (3.7-15)  with  (3.5-5). 

To  evaluate  coefficient  we  use  the  same  procedure  described 
above.  We  omit  details  and  give  only  the  result: 


Htk.  .(-;>wv.. 

ik  (2i+l)i!(2k+l)k!  r=0  s=0 


4I+J+1-S  q=0  : 


r  q/r\  1  /2(j+l-r+q)J 

Z  t*1,l  U/  j+l-r+q"  \  ) 


(3.7-17) 


3.8  Accuracy  Coefficients  for  Rectanqular  Pulses 


If  transmitted  pulses  are  rectangular  with  duration  xp  and  ampli¬ 
tude  p(0),  pulses  of  form  c(t)  will  be  triangular  of  base  duration  2  x 


and  peak  amplitude 

c(0 )  -  p2(0)  x 

r 

according  to  (3.1-3).  Thus, 


(3.8-1) 


c(°)  [1  +  (t/xp)],  -  xp  <  t  5  0 
c(o)  H  -  (t/T  )],  0  <  t  <  tp> 


When  ( 3 . 8-2) i s  substituted  into  (3.3-5)  we  obtain 


(3.8-2) 


3.8-3 
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From  (3.3-7): 


r(-n (2k)! 


(ki )  222k(k+1 ) ! 


Rck=‘ 


Hr  (2k) ! 

(k!)2  22k( k+l)> 


(3.8-4) 


These  functions  are  seen  to  decrease  in  magnitude  rapidly  as  k  in¬ 
creases. 

For  k=0,  R  ^  *  1.  When  (3.8-4)  is  used  in  the  denominator  of  the 
accuracy  coefficient  Rft  given  by  (3.6-2),  this  denominator  behaves  as 
shown  in  Figure  3-2.  as  found  by  computer. 

|  At  this  writing  other  computer  computations  of  the  numerator  of 
R*  and  Rg  [given  in  (3.6-3)]  were  made  but  due  to  slow  convergence  of 
the  series  many  terms  were  found  to  be  necessary  for  accuracy.  Un¬ 
fortunately,  higher  terms  generated  practical  machine  problems  (over¬ 
flows)  that  were  unsolved  as  the  study  ended.  Because  of  these  problems 
further  evaluation  of  R^  and  Rg  remains  for  future  work. 


4.0  ANGLE  ACCURACY  QF 
RHODES'  CLASS  III  MONOPULSE 


4.1  Introduction 

In  a  Rhodes  class  III  system  [2]  having  an  offset  beam  angle 
sensor  (amplitude  monopulse  antenna),  the  sensor  outputs  are  passed 
through  a  short-slot  coupler  to  convert  to  a  phase-sensing  sensor. 

The  coupler  outputs  pass  to  a  receiver  having  limiters  sensitive  only 
to  input  phase  variations.  The  limited  outputs  are  then  compared  in 
an  odd  phase  detector  (one  giving  an  output  proportional  to  the  sine 
of  the  difference  of  the  input  phases)  as  shown  in  Figure  4-1.  The 
odd  phase  detector  can  be  implemented  by  preceeding  the  usual  even 
phase  detector  with  a  90°  phase  shifter  on  one  of  its  inputs. 

A  practical  implementation  of  a  class  III  system,  which  we  shall 
refer  to  as  System  6,  is  shown  in  Figure  4-2  for  a  single  angle 
coordinate. 

In  this  main  section  we  determine  an  expression  for  the  angle 
accuracy  of  the  system  of  Figure  4-2.  We  begin  by  establishing  the 
equivalences  between  outputs  and  P ^  for  a  class  III  system  and 
outputs  £  and  A  of  a  monopulse  system  (1,  p.3]. 

4.2  Amplitude  Monopulse  Equivalences  to  System  6 

An  amplitude  monopulse  antenna  with  its  hybrid  networks  wi’l  pro¬ 
duce  outputs  £,  A  ,  and  A  .  from  its  sum  and  two  difference  channels. 

aZ  ex. 

Within  a  proportionality  constant  that  is  the  same  for  both  monopulse 
and  other  systems,  we  may  consider  Z,  Aa,  and  A  .  as  "patterns"  of  the 

aZ 

antenna. 

It  is  convenient  to  think  of  patterns  £,  A  and  A  „  as  generated 

az  e£. 

by  four  space  patterns  A,  B,  C,  and  D.  Indeed,  when  four  feed  horns 


mixer 


LIM-AMP 


4-1.  Basic  block  diagram  of  Rhodes 
Class  III  monopulse  using  an 
amplitude  monopulse  antenna. 


are  used  with  a  reflector  antenna,  this  is  exactly  the  case  as  shown 
in  Figure  4-3.  Clearly, 
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a2  "  2 


^ 


Pel  '  7  £E  *  ^ 


Pe2  =  2  ‘  AeJ^  * 


(4.2-9) 


(4.2-10) 


(4.2-11) 


When  two  of  these  outputs  (say  and  P ^  corresponding  to  the  azi¬ 
muth  coordinate)  are  processed  using  a  short-slot  coupler  according  to 
Figure  4-5,  the  responses  are  found  to  be 


P,  =  |  [<>/4  E  ♦  e'j*/4  4az] 


(4.2-12) 


P2  ■  \  [eJ’/4  E  -  e-j’/4  4JZ]  • 


(4.2-13) 


I  .  summary,  we  have  found  that  the  responses  P^  and  P 2  for  a 
class  III  system  using  the  angle  sensor  defined  in  Figures  4-4  and 
4-5  are  related  to  the  responses  2  and  A,,  that  the  same  antenna 
would  produce  in  a  sum  and  difference  mode  according  to  (4.2-12)  and 
(4.2-13).  A  similar  statement  can  be  made  for  the  elevation  coordin¬ 
ate  but  following  analysis  is  identical  in  both  coordinates  and  only 
one  analysis  needs  to  be  conducted. 

4.3  Preliminary  Analysis 

In  the  notation  of  Figure  4-2  we  write  (2.2-14)  as 


aZM 


w 

-1  ifi  (R7(t, 

*  loo  Z 

1  /  d  ^{E[Z(t)j} 
\  dx 


{R^t,  t+x)}  dx 


(4.3-1) 


where 


/ 


Fig.  4-5.  Short  slot  coupler  for  use  in 
Rhodes'  Class  III  monopulse. 


EfZ(t)]  =  E[g(t)M(t)]  =  E  [g(t)  V(t)  Y(t)] 


=  E[g(t) 


U(t) 


X(t) 

!x(t)| 


(4.3-2) 


since  the  noises  in  the  two  channels  are  presumed  statistically  in¬ 
dependent.  From  the  analysis  in  Section  3  we  know  the  mean  values  of 
the  limiter  responses  [see  (3.2-6)]  so  (4.3-2)  reduces  to 


E[Z(t)]  »  g(t)  erf 


'  e-j  (t)  ' 

f  e  (t)  - 

erf 

.  /2~  . 

-  /  2  a3  - 

(4.3-3) 
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where  and  a ^  are  the  powers  of  the  noises  n-j(t)  and  n3(t), 
respectively. 

From  the  previous  study  [1]  the  outputs  of  an  amplitude  monopulse 
system  are 


CO 


l  =  1^  G*  l  p(t-xT-nTR)  cos  [uQ(t-TT)  +  eQ] 

n=-°° 

(4.3-5) 

where  K  ,  Gv,  k  ,  etc.  are 
o  l  m 

all  defined  in  [1].  Next,  define 

A  u  r  2 
“i  “  Ko  h 

(4.3-6) 

A  a  p  2 
a2  =  Ko  GS 

k  / 0_  =  a,  k/e, 
m  j  I  m  j 

(4.3-7) 

Per])(t)  = 

CO 

E  p(t  -  XT  -  nTp) 

n=-w 

(4.3-8) 

A  . 

0  =  u  (t  - 

xT)  +  +  (rr/4) 

(4.3-9) 

so 


(1) 

a  . 

Aaz  "  x  a2  Per  (t)  cos  [6  -  (ir/4)] 

(4.3-10) 

*.  \ 

I  *  Per  ^(t)  cos  [0  -  (tt/4)] 

(4.3-11) 

s 

and 

P1  *  2  Per  ^  cos(9)  +  x  sin(0)} 

(4.3-12) 

,  0) 

p2  *  2  Per  ^  {a1  cos(0)  "  a2  x  s1n(9)> 

1 

(4.3-13) 

from  (4.2-12)  and  (4.2-13). 

After  the  matched  filters  In  Figure  4-2  the  signals  become 

e] (0  *  j  Per(t)  s1n(0)  -  «2  x  cos(e)} 

(4.3-14) 

e3(t)  -  ]-  Per(t)  cos(0)  -  a2  x  sln(e)} 

(4.3-15) 

XJ 

where 

CO 

A  — 

•; 

Per(t)  -  1  c(t  -  tt  -  nTp) 

n*  -oo  i  k 

(4.3-15) 

S 

c(t)  *  f  p(0  P  (C  +  t)  d£ 

*  -00 

(4.3-17) 

•- 

c(t)  - ►  |P(u))|2 

(4.3-18) 

4 

P(t)«* - I*  P(w). 

(4.3-19) 

Signals  e^t)  and  e3(t)  given  by  (4.3-14)  and  (4.3-15) 

will  sub- 

sequently  be  used  with  (4.3-3)  in  (4.3-1).  First,  however. 

we  compute 

F*  • 

s 

the  necessary  derivative. 

■» 


€ 
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4.4  Derivative  Calculation 

Since  error  angle  x  does  not  Involve  t  we  can  differentiate  in¬ 


side  the  time  Integral  (time  average) 


d  iA{E[Z(t)]}  « { dE[Z(t)] 


?£]>} 


(4.4-1) 


Lrf  -O/J  .  l_.e 
du  /r 


d  e,(t) 

1 - »  -  (a2/2)  cos(8)  Per(t) 


(4.4-2) 


(4.4-3) 


de3(t> 


and,  if  we  assume 


a  -  (a2/2)  sin(e)  Peru) 


°1  *  °3  *  0 


(4.4-4) 


(4.4-5) 


■  7=ife?=‘fl{9<t)  Per<‘)(s1"(e)  e  63/20  erff/r;] 

*  cos(a)  e'6'  n°  erf  [  !4.< 


+  cos(e)  e  1 


-e  2/2o2 


(4.4-6) 
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After  substituting  the  series  expansions  for  the  error  function  and  the 
exponential  (4.4-6)  can  be  reduced  to 


d  J\  (E[Z(t)]} 

dx  x=0 


00  /  uk  00  ,  .vj,  2k+2£+1 

-/T  a2  r  (-1)  r  (-1)  O] 

110  is  n  k!  (2k+l )  0  _  ~  2k+l .  J~i 

k-0  i"°£!(2/"2  °) 

2k+2£+2  ov+y .  .  2£,.t  2k+2.  2Z 


•  i/Ug(t)[Per(t)]  [si n^K  d(Q)  cos  (0)+cos  (e)  sin  (6)3)  • 

(4.4-7) 

The  factor  involving  sines  and  cosines  varies  rapidly  with  time  com¬ 
pared  with  other  time  factors.  The  time  average  will  then  very  nearly 
equal  the  time  average  of  the  other  factors  multiplied  by  the  constant 
(in  t)  part  of  the  sine/cosine  factor. 

Define 


Cok;  ^  constant  part  of  [sin2k+2(0)  cos 

+  cos2k+2(0)  sin2£(0)] 


(4.4-8) 


and  note  that  in  the  equivalent  amplitude  monopulse  system  the  signal- 


to-noise  ratio  at  the  matched  filter  output  would  be 

Si  ai2c2(0) 


(4.4-9) 


By  using  (4,4-8)  and  (4.4-9)  in  (4.4-7)  we  get 


duA  (E[Z(t)]} 


E  H)*  E  (~^&Cok£ 
7t93  k=0  k!  (2k+l )  £=0  *!22k+2£ 


Per(t)  "1  2k+2£+2'l 


(4.4-10) 
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Two  other  relationships  that  can  be  proved  from  these  two  equations, 
that  are  not  previously  known  to  the  author,  are 


l  (*) 

n=0  ' 


mm)n-  w(\i)£ 


(£+lJ 


(4.6-3) 


mm"-  -k&). 


(4.6-4) 


5.0  SUMMARY  AND  DISCUSSION 


This  study  has  directed  its  efforts  in  four  principal  areas.  The 
first  was  to  rigorously  determine  the  angle  tracking  accuracy  of  an  am¬ 
plitude  comparison  monopulse  radar  having  a  limiter  in  the  sum  channel 
prior  to  the  phase  (angle)  detector.  This  effort  was  successfully 
completed  and  accura:y  is  determined  through  the  variance  of  angle 
error  in  (3.6-1).  The  accuracy  expression  (3.6-1)  requires  computation 
of  accuracy  coefficients  R^  and  Rg,  which  are  given  in  general  by  (3.6-2) 
and  (3.6-3).  These,  in  turn,  require  series  coefficients  and 

'Wrp  35  we"  as  seHes  f’u"ctions  W  Rck  and  RBmHrp:  these  are 
all  given  by  (3.5-2),  (3.5-5)  using  (3.5-6),  (3.5-3),  (3.3-7)  using 
(3.3-5),  and  (3.5-//,  respectively.  All  these  results  are  valid  for  any 
shaped  transmitted  signal,  and  any  duration  range  gate  (range  gate  is, 
however,  limited  to  a  rectangular  shape,  where  the  receiver  is  assumed 
to  have  matched  filters  ). 

The  second  effort  area  was  to  derive  a  general  expression  for  the 
angle  tracking  accuracy  of  any  system,  whether  it  be  linear,  nonlinear 
or  have  any  form  of  signal  processor.  The  main  result  is  (2.2-12), 
which  applies  even  for  off-boresight  tracking  where  xQ  f  0.  For  bore- 
sight  tracking  under  the  special  (but  usually  satisfied)  constraint 
where  the  time-averaged  mean  output  of  the  system  is  an  odd  function  of 
off-boresight  angle  x,  the  applicable  result  is  (2.2-14).  The  result 
(2.2-14)  was  applied  in  the  analysis  described  in  the  foregoing  para¬ 
graph. 

The  third  area  of  effort  centered  around  numerical  evaluation  of 
the  accuracy  coefficients  R^  and  Rg  for  the  amplitude  monopulse  system 
having  a  limiter  and  using  rectangular  transmitted  pulses.  This  effort 


was  only  partially  completed.  Practical  computer  problems  prevented 
evaluation  of  all  but  the  denominator  of  R^,  which  is  plotted  in 
Figure  3-2  for  several  values  of  single  pulse  signal-to-noise  ratio. 

The  fourth  main  area  of  effort  during  the  study  was  to  find  the 
general  accuracy  of  a  Rhodes'  Class  III  monopulse  angle  tracking  sys¬ 
tem.  This  work  was  taken  as  far  as  possible  by  the  end  of  the  study. 
The  interim  result  is  (4.5-3).  Future  work  is  needed  to  obtain  a  final 
variance  expression. 

In  addition  to  the  four  main  work  areas  some  effort  was  made  to 
reduce  the  5-dimensional  sum  required  in  (3.6-3).  Some  progress  was 
made  but  the  work  was  incomplete  at  this  writing.  Other  work  was  al¬ 
so  done  to  model  generally  the  ideal  limiter;  this  work  was  too  incom¬ 
plete  to  include  herein. 
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